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Abstract 


A set of 23,3. r >5 1° x 1° mean free-air anomalies were used to predict 
a set of ff equal area anomalies and their standard errors. Using the 1° 
data incorporating geophysically predicted values of ACIC, 1283 5° blocks 
were computed. Excluding the geophysically predicted anomalies 1249 
blocks were computed. The 1° data was also used to compute covariance 
functions and the equatorial gravity and flattening implied by this data. 

The predicted anomalies were supplemented by model anomalies to form 
a complete 1654 global anomaly field. This data was used in a weighted 
least squares to determine potential coefficients to degree 15, and in a 
summation type formulation to determine potential coefficients to degree 25. 
These potential coefficients sets are compared to recent satellite determinations. 
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1. 0 Introduction 


The combination of satellite and gravimetric data has generally taken 
place using 5°x5° (Rapp, 1967) mean anomalies, or 300 n. m. (Kaula, 1966) 
mean anomalies. In Kaula's work, the 300 n. m. means were estimated on 
the basis of a set of l°xi° mean anomalies provided by the Aeronautical Chart 
and Information Center (ACIC). Using a linear regression technique Kaula 
computed mean anomalies for 935 near equal area blocks. The values for 
934 of these anomalies were given by Gaposchkin and Lambeck (1970). 

An updated set of 1470 5°x5° anomalies was used by Rapp (1969) in 
determining potential coefficients solely from gravity data. The specific 
anomalies were not published although they were used by other investigators 
(e.g. Koch (1970)). 

In defining a new combination procedure (Rapp (1970)) it was necessary 
to determine mean gravity anomalies in 15° equal area blocks. As a first 
step in the evaluation of these anomalies it was decided that a set of 5° or 
300 n. m. blocks should first be computed from the latest available data. 

This paper summarizes the data used, the procedures, the results and the 
analysis of the results, that are related to the determination of these 
anomalies . 

2. 0 The 1° x 1° Free-Air Anomaly Data 


The starting point for these computations was a set of 19, 846 l°xl° mean 
anomalies, referenced to the International Gravity Formula, their accuracy, 
and other pertinent information provided by ACIC (Creighton, 1971). These 
anomalies are a slightly updated set from the 19,164 anomalies described in 
a recent ACIC report (1971). 

The procedures for estimating the l°xi° anomalies were varied but specified 
on the ACIC data. Of the 19, 846 ACIC anomalies 2708 had been computed on 
the basis of gravity-geophysical correlation prediction techniques as described, 
for example, by Durbin (1972). These geophysically derived anomalies are, 
at times, treated separately in this discussion, as they are not based on actual 
gravity measurements. 


Examination of the ACIC material revealed certain areas that were not 
included. It was then decided to do an updating of the ACIC material incorporating 
new material both as additional material and as replacement material. This 
updating process consisted of literature searches and correspondence with 
persons who could make available data that would be helpful to this compilation. 

One of the biggest gaps in the ACIC data occurred in Canada. To obtain 
such anomalies it was necessary to resort to fast, but somewhat less than 
rigorous procedures. In the specific case of Canada, a tape containing 131,105 
point anomalies in Canada was provided to us by J. Tanner. Using this data 
a set of mean Bouguer anomalies were formed for land areas by simple 
meaning the point values that fell within a l°xl° block. These anomalies 
were then converted to free-air anomalies using mean elevations provided 
by ACIC. For ocean areas the mean free air anomalies were computed as 
a mean of the point free-air anomalies. Using this procedure we replaced 323 
of the ACIC anomalies and added 1989 l°xl° anomalies in the Canadian area. 

The procedures followed for other sources depended specifically on the 
source. Some anomalies were obtained directly as recommended values by 
an investigator. Other anomalies were obtained by estimations from anomaly 
contour maps. In many cases we would have more than one anomaly estimate 
for a l°xl° block. The anomaly that was selected to be used was generally 
the one that had the smaller standard error estimated for it. In the case of 
an anomaly given by ACIC and another source, the value given by ACIC was 
accepted if the standard errors were similar. 

At the conclusion of this updating procedure we had examined 11,350 
l°xl° mean anomalies excluding the ACIC anomalies. Of these anomalies, 

862 were used to replace values given by ACIC and 3510 were used where ACIC 
had given no value. Combining this material we obtained a set of 23, 355 l°xi° 
mean anomalies that form the basis for the computations described in this 
paper. 

The accuracy of the l°xl° mean anomalies was determined in several 
ways. The material from ACIC had specific accuracy estimates given. For 
non-ACIC material we automatically assigned a standard error of ±20 mgals 
unless specific information was available on the number of points that were 
used in the estimation of the l°xl° anomaly. If this information was available, 
a standard error was assigned on the basis of the number of points in the 
block. The specific criteria used are given in Table One. 
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Table One 

l°xl° Block Accuracy As a Function of the Number 
of Point Anomalies Used 


No. of Anomalies 


Accuracy 


1-3 
4-7 
8-14 
15 - 25 
26 - 40 
41 - 61 
62 - 85 
86 - 114 
115 - 250 
>250 


± 23 mgals 
21 
19 
17 
15 
13 
11 
9 
7 

± 5 mgals 


The data for Table One was obtained by analyzing the relationship between 
the anomaly accuracy and number of points used in the anomaly estimation 
for all the ACIC anomaly data where such information was given. Such a 
procedure is not ideal because it does not specifically consider the accuracy 
of the observed data nor the distribution of the point data within a block. 

Such considerations would be of prohibitive expense to this study. 


2. 1 Gravity Formula Conversion 


The 1° x 1° anomalies described in section 2.0 were referred to the 
International Gravity Formula. In the computations of this paper we desired 
to refer our anomalies to a gravity formula more consistent with current 
estimates of critical earth constants. Such estimates are (Rapp, 1971): 

a = 6378137. 8 m 

f = 1/298.258 

kM = 3.986013 x 10 14 m 3 /sec 3 

The value of kM includes the mass of the atmosphere so that gravity values 
computed using this kM should (at the surface of the earth) be reduced by 
- 0. 87 mgals (Ecker and Mittemayer, 1967) to obtain a theoretical gravity 
consistent with measurements taken on the surface of the earth, inside the 
atmosphere. The gravity formula consistent with the above constants and 
considerations is: 

y - y t ( 1 + .00530243 sin 3 ® - .00000587 sin 3 2cp) (1) 
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where the equatorial gravity, y E , is 978033.51 mgals. 

In order to convert anomalies given with respect to the International 
Gravity Formula referred to the Potsdam system to the gravity formula 
given in equation (1) which is in an absolute system, we adopt a Potsdam 
correction of -14 mgals (IAG, 1971). If Ag^ are anomalies referred to 
the International Gravity Formula and Agi are the anomalies referred to 
the gravity formula of equation (1), we have: 

Agi = Ag x + ( 1.49 - 13. 71 sin 3 cpy mgals (2) 
where <x> is the mean latitude of the anomaly block. 

Unless otherwise stated, all computations in this paper were done 
using anomalies that had been referred to the gravity formula given in 
equation (1). 

3. 0 Method of Computation of the 300 n. m. Mean Anomalies 


The l°xl° anomalies were processed using basically the same techniques 
and computer programs used by Kaula (1966). The l°xl° means were first 
formed into mean anomalies for areas of 60 n. m. (in latitude) and 60 ±30 n. m. 

(in longitude). These anomalies were then used to estimate the mean anomalies 
in 300 n.m. (in latitude) and 300 ±30 n.m. (in longitude) blocks by predicting, 
in the 300 n.m. block any missing 60 n.m. blocks by linear regression. The 
300n. m. mean anomaly block was then formed as a straight average of the 
25 60 n.m. blocks in the 300 n. m. blocks. No estimations were made for a 
300 n.m. block unless it contained one or more observed 60 n.m. blocks. 

The specific equation used for predicting a 60 n.m. anomaly (g*) was 
given by Moritz (1969): 

g* = S,(C+p)^ £ (3) 

where Cp is a column vector whose elements are the covariance between the 
block (p) to be predicted and the observed anomalies. C is a matrix whose 
elements are the covariances between the observed anomalies; D is an error 
covariance matrix for the known blocks, and g; is a column vector of the observed 
anomalies within the 300 n.m. blocks in which g* was situated. For these 
computations D is taken to be a diagonal matrix with each diagonal element 
being equal to l/m^ where m^ is the accuracy estimate of the observed anomaly 
g 3 . The procedure presented differs from that used by Kaula (1966) in 
that we now do not consider the l°xl° data perfect. 
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The accuracy (m) of the 300 n. m. anomaly was computed from the following 
(Heiskanenand Moritz, 1967, section 7-9, Moritz, 1969, p. 11): 

m 2 = S - C, (C +D)" 1 C t ' (4) 

where C is the mean square value (or variance) of the 300 n. m. mean anomalies, 
and Ci is a column vector representing the covariance between the i th observed 
anomaly and the 300 n. m. block in which it lies. We have: 

n 

Ci=^-£ Ci,* (5) 

k=l 


and 

n n 

(6) 

k = l J=1 

where n is the total number of 60 n. m. blocks within the 300 n. m. block (i. e. n = 25). 

The accuracy estimation obtained through equation (4) will reflect both 
representation errors and errors contributed through inaccurate l°xl° data. 


4.0 300 n. m. Anomaly Set 


4. 1 Numerical Procedures and Results 

We first determined the covariances from the 60 n. m. block anomalies 
that are needed in the evaluation of the C lf k and Cj values of equations 
(5) and (6). For comparison purposes this was done using all the 1° x 1° data 
and then a 1° x 1° data set that excluded the geophysically predicted anomalies. 
The covariances were computed using the same procedures and program used 
by Kaula (1966). The results are given in Table Two. 
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Table Two 

Short-Range Autocovariances of 60 ± 30 n. m. 
Free Air Gravity Anomalies 


Geophysically Predicted 

All Anomalies Used Anomalies Excluded 


stance 

Number of Pairs 

Covariance 

Number of Pairs 

Covariance 

deg 

in Sample 

mgal 8 

in Sample 

mgal 2 

0.00 

19837 

690 

19245 

697 

0.97 

22592 

317 

21198 

309 

1.60 

41526 

227 

38956 

217 

2.50 

55997 

163 

52182 

147 

3.46 

42947 

95 

40318 

71 

4.37 

24849 

39 

23167 

17 

5.34 

3810 

-16* 

3579 

-42* 

6.18 

110 

-63* 

110 

-92* 


* interpolated values 

, Using the 1° covariances where all the data was used a set of covariances 
Ci and C_i , C were computed and used to carry out the prediction of 1283 
300 n. m. blocks based on all the 1° x 1° data._The anomaly variance from 
the predicted anomalies was 245 mgal 2 while C was 182 mgal 2 as computed 
from the 1° covariance data. This latter value is smaller than that actually 
found due to the fact that more anomaly variability exists in the actual 
anomalies than is represented by the covariance functions. In order to 
use the more realistic higher values of Ci , Cj, and C the covariance 
functions previously obtained for these quantities were scaled by a factor 
which would make C equal to the variance of the 300 n. m. means as obtained 
from the first prediction tests. Using these covariances, the 1283 anomalies 
were predicted again, as well as a set of 1249 300 n.m. anomalies that 
wrer predicted excluding all geophysically predicted 1° anomalies. Maps 
showing the location of the 1283 anomaly set are given in Figure 1A and IB 
where the location of the anomaly is indicated by the integer part of the rounded 
standard error of the prediction. 

The results of these predictions are given in Table Three. The first 
section of the table given the 1249 anomalies predicted excluding all geophysically 
predicted anomalies. These anomalies are followed by 121 anomalies which 
were predicted from the 1° anomaly set that included the geophysically predicted 
anomalies. Thirty-four of the 121 values do not appear in the 1249 anomaly 
set. If one wishes to have a set of 300 n.m. mean anomalies based only on 
actual gravity observations, the 1249 set should be used. If a requirement 
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Figure 1A 

Standard Errors of 1283 Predicted Anomalies 
Northern Hemisphere 
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Figure IB 

Standard Errors of 1283 Predicted Anomalies 
Southern Hemisphere 
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Table Three 

1249 300 n.m. free-air mean anomalies computed using 
no geophysically predicted 1° x 1° anomalies 
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exists for data that incorporates the geophysically predicted anomalies, 87 
of the additional 121 anomalies should be used as replacements in the 1249 
set while the remaining 34 new anomalies should be added to the 1249 set. 

The specific data in the table consists of the latitude and longitude 
of the center of the block, the anomaly (mgals), the anomaly standard 
error, and the number (n) of known 60 n. m. blocks used in the computation 
of the 300 n.m. block. All anomalies and their standard errors have been 
rounded to the nearest mgal. Consequently there are several blocks that 
have a standard error of zero printed (and given in Figure 1) since the 
actual standard error was less than 0. 5 mgals. These blocks and the actual 
standard errors are: co = 47. 5?, A=11.0°, m = ± 0. 4 mgals; cp = 42. 5°, 

A = 268. 5°, m = ±0.4 mgals; 9 = 37. 5°, A= 275.0°, m =±0.3 mgals. 


4. 2 Comparison with Kaula (1966) Anomalies 


The values of 934 anomalies used by Kaula (1966) are given in Gaposchkin 
and Lambeck (1970). These anomalies were compared to the corresponding 
anomalies of the 1283 set for various ranges of accuracy of the anomalies of 
the 1283 set. The results of these comparisons are shown in Table Four. 


Table Four 

Comparison of Kaula Anomalies with 
Anomalies from 1283 Set 


Maximum 

Number of 

Mean Difference 

Root Mean Square 

Accuracy (mgals) 

Blocks 

(1283 Set - Kaula) 

Difference 

1 

3 

-1. 5 mgals 

±3.8 mgals 

2 

37 

-3.7 

5.6 

3 

116 

-4.4 

8.2 

4 

212 

-2.9 

12.3 

10 

697 

-2.3 

11.5 

16 

919 

-2.0 

10.9 


We can see that there is a systematic anomaly difference between the 
Kaula results and the results of this paper. This difference probably represents 
the use of different reference equatorial gravity values. 

These comparisons also revealed that the Kaula anomaly set contained 
15 anomaly values that were not contained in the 1283 set. These values are 
given in Table Five. 
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Table Five 

Kaula Anomalies Missing from 1283 Set 


lO 

$ 

X° 

Ag (mgals) 

n 

m (mgals) 

82. 5 

100.0 

-4 

1 

± 14. 5 

67. 5 

70. 5 

-19 

1 

14.5 

57. 5 

337.0 

-9 

5 

11.1 

52. 5 

143.0 

-3 

4 

11.8 

22.5 

169. 5 

-8 

5 

11.1 

12.5 

249. 5 

7 

8 

9.1 

7. 5 

164.5 

4 

4 

11.8 

7. 5 

337. 5 

5 

8 

9.1 

2.5 

297. 5 

7 

1 

14. 5 

2. 5 

337. 5 

-2 

8 

9.1 

-7. 5 

210.5 

14 

1 

14.5 

-12. 5 

208. 5 

16 

6 

10.4 

-17. 5 

2. 5 

6 

3 

12.7 

-17. 5 

352. 5 

42 

2 

13.2 

-87. 5 

60.0 

-17 

5 

11.1 


In this table n is the number of 60 n. m. blocks used in the prediction and 
m is the standard error of the prediction based on n as determined from 
corresponding predictions made in the 1283 set. 

These values could be incorporated into our 1283 anomaly set by taking 
into account the systematic difference indicated in Table Four. 


4. 3 Long Range Covariances 


The 300 n. m. predicted mean anomalies (with and without geophysically 
predicted data) were used to determine covariance functions for this block 
size. Using programs originally written by Kaula, these covariances were 
computed and are given in Table Six. These values may be compared to 
similiar values given in Table Two of Kaula (1966). 
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Table Six 

Long-Range Autocovariances 
of 300 n. m. Free-Air Gravity Anomalies 
(mgal 3 ) 


Distance 

All Anom. 

Geophy. Anom. 

Distance 

All Anom. 

Geophy. Anom. 

(deg) 

Used 

Excluded 

Cdeg) 

Used 

Excluded 

0 

250 

245 

90 

-7 

-8 

5 

139 

133 

95 

-2 

-2 

9 

102 

99 

100 

2 

3 

13 

63 

63 

105 

7 

8 

18 

39 

40 

111 

10 

12 

23 

21 

20 

116 

11 

14 

29 

6 

5 

121 

11 

14 

34 

-2 

-1 

126 

9 

10 

39 

-4 

-4 

131 

6 

7 

44 

-5 

-4 

136 

3 

2 

49 

-8 

-7 

141 

0 

-2 

54 

-10 

-11 

146 

-3 

-3 

59 

-12 

-14 

151 

-2 

-3 

64 

-9 

-14 

156 

-1 

-1 

69 

-9 

-12 

162 

-3 

-4 

74 

-8 

-10 

167 

-10 

-10 

80 

-7 

-10 

172 

-18 

-13 

85 

-8 

-10 

175 

-20 

-17 


4. 4 Anomaly Degree Variances from Covariance Data 

The anomaly degree variances can be computed from covariance data 
from the following well known expression: 

TT 

0 "f(Ag) = — J ( cost ) C( ^ ) sin d<|r (7) 

In this equation ijr is the spherical arc separation, P ^ is a Legendre polynomial 
of degree & , and C(i(r) is the covariance function. The results of applying 
equation (7) to the covariance values given in Table Six are given in Table Seven 
along with values from the gravimetric data as given by Kaula (1966, Table 
Three) and values obtained from satellite derived potential coefficients as 
given by Gaposchkin and Lambeck (1970). 
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Table Seven 



Anomaly Degree Variances (^(Ag) 
(mgal 2 ) 

Gravimetric 

Satellite 



C(t) Table 5 

C(i|0 Table 5 


l 

Kaula (1966) 

All Data 

No Geophy. 

Gaposchkin & Lambeck (1970) 

0 

2.7 

1.2 

0.8 

— 

1 

-0.5 

2.2 

1.2 

— 

2 

6.3 

12.7 

13. 1 

7.4 

3 

31.8 

31.3 

34.0 

33.3 

4 

18.6 . 

13.6 

13.0 

19.7 

5 

8.4 

15.1 

11.4 

17.5 

6 

22.2 

19.9 

20.4 

14.4 

7 

11.0 

15. 5 

14.8 

16.4 

8 

9.2 

7.5 

8.3 

8.5 

9 

10.1 

16.1 

16.3 

15.1 

10 


9.6 

9.1 

17.7 

11 


10.8 

8.9 

13.7 

12 


3.9 

3.5 

8.4 

13 


8.2 

6.2 


14 


8.2 

8.7 


15 


8.1 

.8.7 



These values will be compared in a subsequent section to anomaly degree 
variances computed from potential coefficient solutions made using the gravity 
material computed in this work. 


5. 0 Gravity Formula Parameters 


The 1° x 1° gravity data and the covariance functions previously computed 
enable us to determine values of equatorial gravity and the flattening as implied 
by terrestrial gravity material. To determine both equatorial gravity and the 
flattening we perform a least squares adjustment that minimizes the sum of 
the squares of the weighted anomalies using as a model the gravity formula 
in the form of equation (1). Such a computation was carried out with the 20, 662 
1° x 1° anomalies that did not contain any geophysical anomalies using several 
weighting schemes and with a constraint on the flattening forcing it to be 
equal to 1/298.258. These solutions are described as follows: 

Solution One: y £ and f adjusted, weights = coscp/m 2 
Solution Two: f constrained, weights = cos cp / m 2 
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Solution Three: y E and f adjusted, weights = cos 9/ 10 3 
Solution Four: f constrained, weights = cos 9/ 10 3 

Solutions Three and Four assume the anomaly data has the same accuracy. 

In practice this reduces the dominance of the more accurately observed 
land areas. 

In addition to using the gravity formula technique for finding y E , we 
may use the zeroth order anomaly degree variance given in Table Seven as 
computed from the covariance function where no geophysically predicted 
anomalies were used. We have (Heiskanen and Moritz, 1967): 

/ 1 

v Vo (Ag) = Ago = ( y s - y«) - — y( fu - u ) (S) 

3 

where y 8 , f 0 are equatorial gravity and flattening of the "best" parameters 
while y„, f« are corresponding parameters of the gravity formula to which 
the anomalies have been referred. Since the anomalies, and thus the anomaly 
covariances and the anomaly degree variances have been referred to a gravity 
formula where f is considered to be accurately known, we have ^ = f R . 

Thus : 


y& = Tr +V (To (Ag) (9) 


This computation of equatorial gravity constitutes Solution Five. 

Solution Six for y E and f are the values and constants inherent with equation (1). 
Since these constants have been derived independently of gravity information, 
their comparison with the gravity dependent results is of interest. 

The results of these solutions are summarized in Table Eight where all 
equatorial gravity values have been referenced to an absolute system using 
a Potsdam correction of - 14 mgals. 
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Table Eight 

Equatorial Gravity and Flattening Determinations 


Solution 

y E (mgals) 

1/f 

One 

978034.4 

298.646 

Two 

978035.8 

298.258* 

Three 

978032. 1 

298.265 

Four 

978032.1 

298.258* 

Five 

978034.4 

298.258* 

Six 

978033.5 

298.258* 


Comparing solutions of the same flattening, we find solution 5, followed by 
solution 4 agrees best with solution 6. Solutions 2 differs by 2. 3 mgal from 
solution 6 indicating some distortion has taken place by weighting according 
to actual anomaly accuracy estimates. We could conclude that within the 
accuracy of the data there is reasonable agreement between equatorial gravity 
values determined from gravity data and satellite data. 


6.0 Potential Coefficient Determinations 


The 300 n. m. mean anomalies computed in this paper may be used to 
determine potential coefficients describing the earth’s gravitational field. In 
order to do this we first consider the relationship between point gravity anomalies, 
referred to a given gravity formula and potential coefficients: 


Agr = Ag t 


kM 


I (C i. c03mX+s *. sin mX) P^jsincp) J 

( 10 ) 


4=2 


8=0 


where Ag 0 is the mean global value of Ag r , r is the geocentric distance to the 
computation point, C* and S^ a are fully normalized potential coefficients where 
C* represents the difference between the actual values and reference values 
consistent with the flattening to which the anomalies are referrenced. P is 
the fully normalized Legendre functions and 9 is the geocentric latitude. 


Generally equation (10) has been evaluated using center point block 
coordinates. With higher degree solutions being continually made, errors will 
start to exist at the higher degrees unless an integrated form of equation (10) 
is used that recognizes that the anomaly values are given in blocks bordered 
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by meridians and parallels, e. g. cpi, <ps, X lf X 2 . Thus a more accurate expres- 
sion relating potential coefficients to mean anomalies is as follows (Rapp, 1972, 
Desrochers, 1971): 


Ag- = Ag 0 + 


Ts 


kM 


AX(sincp 2 - sincpx ) f 


- e 


[i ii - i > (t)*[ Axe *^“' 


1 = 2 


*J P^ j0 (sincp) cos cp dcp + B (sinmX a - sinmXj ) -S^cos mX a - 


= i 




cosmXx) ] J P^ a (sincp) coscp dcp J J (11) 

r to. 


where r is the average geocentric radius in the block which, in these computa- 
tions, was taken at the mean geocentric latitude of the block, and k is the highest 
degree to which the summation is to be carried, or for which a potential coefficient 
solution is to be made. 

Although the complications resulting in going from equation (10 ) to equation 
(11) are not justified for the lower degree solutions, it is a necessity as solutions 
are generated for the higher degrees, say 15 or more. 

In solving for potential coefficients from gravity material it is desirable 
to have a global set of anomalies which, in the Kaula divisions used in this paper, 
will consist of 1654 blocks. To obtain such a global set the model anomalies 
of Uotila (1964) were used to supplement the predicted anomaly sets assigning 
such model anomalies a standard error of ± 20 mgals as has previously been 
done (Rapp, 1969 ). 

Using equation (11) as a model several weighted least squares solutions 
were made. In applying equation (11 ) , the value of r was computed using the 
value of a = 6378137.9 m and f = 1/298. 258 and all block division latitudes 
defined in terms of geodetic latitudes were converted to geocentric latitudes. 

The longitude limits on the blocks were individually computed for each block 
based on the specific 1 anomaly selection criteria of the anomaly prediction 
program. 

In addition to the relationships represented by equation ( 10 ) (or (11)), 
potential coefficients may be directly computed from the anomaly data from 
the following expressions: 
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i r r a - r cos mX i 

I^TTT n **l'C*« UmX W < 12 > 

This formula may also be written in a manner that obtains a mean value of the 
kemal (excluding Ag) in a manner analogous to what was done in writing 
equation (11). If, in equation (10) (or (11)), r was set equal to a, kM/r 3 was 
set to y, and a least squares adjustment was made with equal anomaly variances, 
the results for the potential coefficients from either (11) or (12) would be the 
same. Consequently, the application of equation (12) may be viewed as a 
least squares adjustment with equal anomaly variances. The computation of 
potential coefficients through equation (12) is considerably faster than the 
procedure carried out using equation (10) or (11) because in the latter case, 
normal equations and their solutions need to be carried out. 

6. 1 Solutions Made 

Several potential coefficient solutions were made with the following 
being of primary interest: 

Solution One: Complete to degree 12 using 1283 anomalies based on 

predictions incorporating the 1° geophysically predicted 
data, plus 371 model anomalies 

Solution Two: Complete to degree 12 using 1249 anomalies based on 

predictions excluding the 1° geophysically predicted 
data, plus 405 model anomalies 

Solution Three: Complete to degree 12 using the 934 Kaula anomalies 
plus 720 model anomalies 

Solution Four: Complete to degree 12 as given by Rapp (1969) 

Solution Five: Complete to degree 15 using same data as solution one 

Solution Six: Complete to degree 25 using same data as solution one 

but using equation (12) for the computation 

Solution Seven: Complete to 15,15 using 1283 anomalies plus the 15 
Kaula anomalies (Table 6) missing in the 1283 set, 
plus 356 model anomalies 


j \ 
K j 
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In Solutions One, Two, Three, and Five the value of Ag 0 was computed as 
the mean value of the 1654 anomaly set. The potential coefficients and their 
standard errors from Solution Five, and the potential coefficients from Solution 
Six are given in Table Nine. The geoid undulations implied by the potential 
coefficients of Solution Five with respect to f = 1/298. 258 and a best fitting 
"a" are given in Figure Two as computed by Method One in Rapp (1971). 

6. 2 Comparison of Solutions to Satellite Derived Potential Coefficients 


In order to judge the accuracy of these solutions, we compare them to 
the recent satellite derived potential coefficients of Gaposchkin (1969) and 
Smith, Lerch, and Wagner (1972). The comparisons are made with respect 
to three quantities: 1. the root mean square coefficient difference; 2. the 
solution correlation (r) as defined by Rapp (1972); 3. the percentage difference 
(p). This latter quantity is computed as the average percentage difference 
of each degree, where the percentage difference of each degree (p ^) is computed 
from the following expression : 

p. = 100 • J — -i ZA3 (13) 

1 V 24+ 1 £cy 3 ’ 


where a is the value of a potential coefficient in degree l of the satellite solution. 
These results are given in Table Ten. 


Table Ten 

Comparison of Potential Coefficient Sets 


Comparison Coefficients 


Solution 

Gaposchkin (1969), 
to 8 only 

A x 10 s r p(%) 

Smith et als. (1972), 
to 8 only 
A x 10 s r p(%) 

Smith et als. (1972) 
to 12 only 
A x 10 6 r p(%) 

One 

.190 

.937 

16.0 

. 187 

.940 

14.4 

.140 

.928 

18.6 

Two 

.191 

.934 

15. 7 

. 190 

.936 

14.5 

.141 

.925 

18.3 

Three 

.344 

. 762 

24. 1 

. 339 

.790 

21.8 

.236 

. 778 

22.2 

Four 

.239 

.895 

16.4 

.246 

.892 

15.9 

. 173 

.886 

17.8 

Five 

. 188 

.938 

15.3 

. 185 

.942 

13.9 

.134 

.934 

16.3 

Six 

. 186 

.938 

14.5 

. 180 

.947 

12.9 

.128 

.942 

14.4 

Seven 

.192 

.937 

15.6 

. 188 

.940 

14.1 

.136 

.933 

16.4 
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Table Nine 

Potential Coefficients From Gravity Data (X 10 s ) 


Solution Five Solution Six 


l 

m 

C £a 

m (C) 


m(S) 



? 

0 

-484.003 

0.285 



-484.400 


2 

2 

3.177 

0.327 

-1.635 

0.308 

2. 721 

-0.828 

3 

0 

0. 749 

0.148 



0.453 


3 

1 

L. 290 

0.154 

-0. 189 

0.134 

1. 224 

0.069 


2 

1.020 

0.165 

-0.603 

0.163 

0.998 

-0.384 

3 

3 

0.664 

0. 161 

1.305 

0.169 

0.794 

1.159 

4 

0 

0.662 

0.105 



0.614 


4 

l 

-0.288 

0.111 

-0.210 

0.105 

-0.348 

-0.217 

4 

2 

0.339 

0. 115 

0. 244 

0. 123 

0.413 

0.281 

4 

3 

0.336 

0.104 

-0.401 

0.109 

0. 758 

-0.198 

4 

4 

-0.053 

0.112 

0.268 

0.112 

-0. 084 

0.237 

5 

0 

0.001 

0.083 



-0.078 


c, 

J 

1 

-0.105 

0.082 

-0.027 

0.074 

-0. 192 

-0.133 

5 

2 

0.521 

0.086 

-0.085 

0.090 

0.458 

-0.067 

5 

3 

-0.350 

0.086 

-0.003 

0.090 

-0. 245 

-0.143 

5 

4 

0.086 

0.079 

0.003 

0.079 

0.025 

-0.011 

5 

5. 

0.286 

0.084 

-0.681 

0.087 

0. 120 

-0.461 

6 

0 

0.087 

0.065 



-0.003 


6 

l 

-0.068 

0.067 

-0.026 

0.063 

-0.059 

-0.039 

6 

2 

0.193 

0.065 

-0. 178 

0.066 

0. 149 

-0.170 

6 

3 

-0. 099 

0.071 

-0.110 

0.075 

-0. 109 

-0.052 

6 

4 

-0.205 

0.069 

-0.419 

0.068 

-0. 130 

-0.332 

6 

5 

-0.394 

0.062 

-0.591 

0.066 

-0.311 

-0.451 

6 

6 

0.007 

0.070 

-0.230 

0.068 

-0.013 

-0.137 

7 

0 

0.114 

0.054 



0.092 


7 

1 

0. 159 

0.057 

0.180 

0.052 

0.228 

0.145 

7 

2 

0.291 

0.052 

0.058 

0.054 

0.273 

0.077 

7 

3 

0.099 

0.056 

-0. 123 

0.060 

0. 134 

-0.096 

7 

4 

-0.158 

0.061 

-0.13? 

0.060 

-0. 143 

-0.127 

7 

5 

-0.035 

0.056 

0. 048 

0.057 

-0.036 

0.055 

7 

6 

-0.240 

0.053 

0.157 

0.053 

-0. 189 

0.135 

7 

7 

-0.049 

0.058 

-0.109 

0.057 

-0.011 

-0.031 

P 

0 

-0.045 

0.046 



-0.004 


8 

l 

-0.137 

0.048 

0.070 

0.044 

-0.056 

0.039 

8 

2 

0.001 

0.046 

0. 148 

0.047 

0.054 

0.164 

8 

3 

0.083 

0.046 

0.008 

0.047 

0. 059 

0.015 

8 

4 

-0.160 

0.053 

-0.002 

0.051 

-0.099 

0.039 

8 

5 

0.009 

0.050 

0.043 

0.052 

0.002 

0.041 

8 

6 

-0.139 

0.047 

0. 130 

0.046 

-0.070 

0.122 

8 

7 

0.018 

0.047 

0.094 

0.045 

0.041 

0.120 

8 

8 

-0.103 

0.050 

0.047 

0.049 

-0.082 

0.039 
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9 

0 

0.174 

0.040 



0 . 118 


9 

l 

0.068 

0.042 

0.019 

0.039 

0 . 142 

- 0.003 

9 


0.137 

0.040 

- 0.145 

0.041 

0.061 

- 0.075 

9 


- 0.205 

0.040 

- 0.082 

0.041 

- 0 . 162 

- 0.042 

CJ 

4 

- 0.072 

0.043 

- 0.020 

0.041 

- 0. 043 

0.023 

9 

■5 

- 0.057 

0.045 

0.005 

0.047 

- 0.094 

0.043 

9 

6 

0.072 

0.043 

0 . 175 

0.043 

0.034 

0.147 

9 

7 

- 0.017 

0.040 

0.078 

0.040 

- 0.022 

0.028 

9 

8 

0.262 

0.041 

0.023 

0.040 

0 . 201 

0.024 

9 

9 

0.052 

0.043 

0.036 

0.044 

- 0.005 

0.035 

10 

0 

0.020 

0.035 



0.011 


10 

1 

0.047 

0 . 038 

- 0.023 

0 . 03 5 

0.054 

- 0.056 

10 

'2 

- 0.057 

0.035 

- 0.122 

0.036 

- 0.057 

- 0.048 

10 

3 

- 0.014 

0.035 

- 0.082 

0.037 

- 0.028 

- 0.069 

10 

4 

- 0.070 

0.037 

- 0.096 

0 . 036 

- 0.063 

- 0.102 

1 0 

5 

- 0.048 

0.038 

0.014 

0.039 

- 0.032 

- 0.007 

10 

6 

- 0.006 

0.040 

- 0.094 

0.040 

- 0.011 

- 0.055 

10 

7 

0.135 

0.038 

- 0.070 

0.037 

0.082 

- 0.013 

10 

8 

0.006 

0.036 

- 0.155 

0.035 

0.012 

- 0.097 

10 

9 

0.096 

0.036 

- 0.010 

0.036 

0 . 133 

0.004 

10 

1.0 

0 . 150 

0.038 

- 0.002 

0.039 

0.047 

0.001 

1 1 

0 

- 0.076 

0.031 



- 0.066 


11 

1 

- 0.027 

0.03 3 

0.046 

0.030 

- 0.005 

0.017 

11 

2 

- 0.011 

0.032 

- 0.131 

0.033 

0.000 

- 0.082 

11 

3 

- 0.078 

0.031 

- 0.072 

0.032 

- 0.075 

- 0.040 

li 

4 

- 0.082 

0.033 

- 0 . 133 

0.0 32 

- 0.060 

- 0.091 

1 1 

5 

0.061 

0.033 

0.034 

0.034 

0.063 

0.009 

11 

6 

- 0.012 

0.037 

- 0.071 

0.035 

0.039 

- 0.045 

11 

7 

0.042 

0.036 

- 0.122 

0.035 

0.016 

- 0.044 

11 

9 

- 0.058 

0.033 

0.084 

0.032 

- 0.022 

0.062 

1 1 

9 

0.011 

0.03 1 

0.030 

0.031 

0 . 009 

0.060 

11 

10 

- 0.065 

0.032 

- 0.065 

0 . 032 

- 0.042 

- 0.046 

11 

11 

0.136 

0.035 

- 0.032 

0.035 

0.055 

0.020 

1 2 

0 

- 0.042 

0.029 



- 0.045 


1 ? 

1 

- 0.067 

0.029 

0.012 

0.026 

- 0.030 

- 0.043 

1 ? 

2 

- 0.067 

0.029 

0.038 

0 . 030 

- 0.038 

0.007 

1 ? 

3 

0.071 

0.028 

0.027 

0. 029 

0 . 026 

- 0.018 

1 ? 

4 

- 0.021 

0.02 9 

- 0.020 

0.028 

- 0.046 

- 0.023 

i ?. 

5 

0.052 

0.029 

- 0.008 

0.029 

0.040 

0.007 

12 

6 

0.001 

0.031 

0.054 

0.030 

0.032 

0.037 

12 

7 

- 0.079 

0.032 

0.023 

0.031 

- 0.052 

0.017 

12 

3 

0.015 

0.030 

0 . 059 

0.030 

0 . 003 

0.022 

12 

9 

- 0.056 

9.02 8 

0.049 

0.028 

- 0.021 

0.013 

12 

10 

0.004 

0.02 8 

- 0.058 

0.027 

- 0.008 

0.001 

12 

t 1 

0.024 

0.02 9 

- 0.008 

0.029 

- 0.009 

- 0.004 

12 

12 

0.020 

0.032 

- 0.038 

0.031 

0.031 

- 0.001 
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1 ? 

0 

0.031 

0.025 



0.036 


13 

l 

0.024 

0.026 

- 0.037 

0.023 

0.009 

- 0.030 

13 

2 

- 0.032 

0.025 

- 0. 056 

0.025 

- 0.022 

- 0.025 

13 

3 

- 0.008 

0.026 

0.031 

0.026 

0 . Oil 

0.036 

13 

4 

0.029 

0.026 

- 0. 073 

0.025 

- 0. 007 

- 0.030 

13 

5 

0.071 

0.026 

0 . 07 ? 

0.026 

0.070 

0.062 

1.3 

6 

- 0.065 

0.028 

0.080 

0.026 

- 0.061 

0.011 

13 

7 

- 0.093 

0.028 

0.071 

0.028 

- 0.031 

0.034 

13 

8 

0.003 

0.028 

0.04 0 

0.02 8 

- 0. 009 

- 0.005 

1.3 

9 

0.003 

0.026 

0.045 

0.027 

- 0.002 

0.058 

13 

10 

0.020 

0.024 

0. 029 

0.024 

- 0. 008 

0.012 

13 

11 

- 0 . 0 R 5 

0.024 

- 0.054 

0.024 

- 0.013 

0.006 

13 

12 

- 0. 032 

0.026 

0.081 

0.026 

0.006 

0.073 

13 

13 

- 0.089 

0.028 

0.125 

0.028 

- 0. 056 

0.048 

14 

0 

0.058 

0.022 



0 . 014 


14 

1 

- 0.046 

0.024 

0 . 035 

0.021 

- 0.002 

0.023 

14 

2 

- 0.020 

0.022 

0.040 

0.022 

- 0.024 

0.009 

14 

3 

0.039 

0.024 

- 0.031 

0.02 3 

0.022 

- 0.022 

14 

4 

- 0.047 

0.024 

0.022 

0.023 

0.015 

0.014 

14 

5 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

14 

6 

- 0.015 

0.025 

- 0.056 

0.024 

0.010 

- 0.003 

14 

7 

0.062 

0.025 

0.014 

0.024 

0 . 022 

- 0.025 

14 

8 

0.013 

0.025 

- 0.060 

0.025 

- 0.029 

- 0.033 

14 

q 

0.020 

0.025 

0.025 

0.024 

- 0.007 

0.075 

14 

10 

0.035 

0.022 

- 0.074 

0.022 

0.048 

- 0.036 

14 

11 

0.046 

0.021 

- 0.068 

0.021 

0.017 

- 0.037 

14 

12 

- 0.023 

0.021 

- 0.026 

0.021 

0.009 

- 0.036 

1.4 

13 

0.0 53 

0.022 

0.016 

0.023 

0.001 

0.069 

14 

14 

- 0.005 

0.024 

- 0.001 

0.024 

- 0.035 

- 0.014 

1 5 

0 

0.013 

0.018 



0 . 016 


15 

l 

0.001 

0.019 

0.079 

0.017 

0.039 

0.015 

15 

7 

- 0.050 

0.018 

- 0.038 

0.018 

- 0. 019 

- 0.000 

15 

3 

- 0.012 

0.018 

0.084 

0.019 

0.012 

0.038 

15 

4 

0.016 

0 . 018 

0.021 

0.018 

- 0.013 

- 0.010 

15 

5 

- 0.025 

0.018 

0.046 

0.017 

0.026 

0.008 

1 5 

6 

- 0.008 

0.018 

- 0.120 

0.018 

0.026 

- 0.065 

1 5 

7 

0.061 

0.018 

0 . 049 

0.018 

0.045 

0.027 

l 5 

a 

- 0.076 

0.018 

- 0.017 

0.018 

- 0. 024 

0.012 

15 

9 

- 0.070 

0.018 

0.047 

0.018 

0.002 

0.030 

15 

10 

- 0.001 

0.018 

0.028 

0.018 

- 0.017 

- 0.003 

15 

11 

- 0.050 

0.017 

- 0.005 

0.018 

0.015 

0.005 

15 

17 

0.012 

0.017 

0.043 

0.017 

- 0.012 

0.946 

15 

1 3 

- 0.005 

0.018 

0.000 

0.018 

- 0 . 041 

- 0.004 

15 

1 4 

- 0.027 

0.019 

- 0.046 

0.019 

- 0.002 

- 0 . 01 7 

15 

15 

- 0.092 

0.020 

0.046 

0.021 

- 0.041 

0.028 
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From the table we see the following: 


1. The potential coefficient solutions based on gravity data agree slightly better 
with the new satellite derived solution of Smith et als. than with the Gaposchkin 
solution. 

2. The poorest solution is Solution Three where the Kaula anomalies were used. 
The second poorest solution is the older solution of Rapp (Solution Four). 

The better agreement of the newer solutions with the satellite derived 
coefficients reflects the increase and improvement of the current gravity 
material over that used in the earlier solutions. 

3. There is very little difference between the coefficient sets that include 
and exclude the geophysically computed 1° anomalies in their formation 
(Solutions One and Two). 

4. The inclusion of the Kaula derived anomalies into the 1283 set (Solution 

Seven) appears to have slightly deteriorated the potential coefficient sets 
when compared to the results of Solution Six. 

5. Of the weighted least squares solutions made, the best set of potential 
coefficients appears to be those found in Solution Five. 

6. Of all solutions, the best agreement with the satellite derived solutions is 
with the coefficients of Solution Six. Information to be considered in the 
next section, however, reveals some undesirable features of Solution Six. 
Consequently for the results of this paper, the potential coefficients from 
Solution Five will be considered the best overall set. 


6. 3 Anomaly Degree Variances from Potential Coefficients 


The anomaly degree variances discussed in section 4.4 may also be 
computed from potential coefficients using the following equation: 

Z 

o*(Ag) - y 3 (Z- 1) 2 £ (C- + s; b ) (12) 

m — 0 

where y is a mean value of gravity. Such values have been computed from several 
solutions and are given in Table Eleven along with the anomaly degree variances 
given in Table Seven as computed from the covariance data obtained from the 1283 
anomaly set. 
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Table Eleven 

Anomaly Degree Variances from Potential Coefficient Solutions 

(rngal 2 ) 


Solutions 

From Table 6 
(Col. 3) 

l 

One 

Five 

Six 

2 

11.3 

12.3 

7.8 

12.7 

3 

24.5 

22.3 

18.5 

31.3 

4 

11.6 

10.8 

9.7 

13.6 

5 

13.2 

14.8 

9.0 

15.1 

6 

22. 7 

21.1 

12.4 

19.9 

7 

16.6 

11.3 

9.9 

15.5 

8 

9.1 

6.6 

4.5 

7.5 

9 

16.2 

15.3 

9.3 

16.1 

10 

14.3 

10.3 

5.6 

9.6 

11 

14.9 

12.3 

5.4 

10.8 

12 

15.7 

5.3 

2.2 

3.9 

13 


13.0 

4.9 

8.2 

14 


7.9 

4.1 

8.2 

15 


12.9 

4.0 

8.1 


The larger degree variances for the higher degrees for Solution One 
reflect the aliasing effect 'Desrochers, 1971) caused by the truncation of the 
of the least squares adjustment. This truncation causes a distortion in the 
coefficients of the degrees near the truncation degree. This effect is also 
noticable in Solution Five, but to a lesser extent. The degree variances 
from Solution Six are all smaller than the ones from Solution Five and those 
computed from the covariance function data. The unreasonably small values 
of the anomaly degree variances from Solution Six indicate that solution to 
be undesirable as compared with other solutions (such as Five). 


7. 0 Summary 


This paper had described the estimation of a 300 n. m. (or 5° equal 
area) block terrestrial free air anomaly field. The derived fields were 
analyzed to obtain covariance functions, anomaly degree variances, and 
potential coefficients. Only small and insignificant differences between 
solutions that did and did not incorporate the geophysically estimated 1° x 1° 
into the 300 n. m. predicted set were found. 

The potential coefficients as determined from this new data agree ' 
better with the satellite derived potential coefficients than the coefficients 
of earlier sets derived from gravity data. This improvement reflects the 
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improvement in the gravity coverage and the accuracy of the material since 
the earlier solutions. 

The gravity data of this paper will be used for the formation of Iff and 
perhaps 10° equal area block mean anomalies. In addition a combination of 
this gravimetric data with satellite derived potential coefficients will be 
carried out to determine a 300 n.m. gravity field that is consistent with 
a high degree potential coefficient field. 
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